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1. Introduction 

The observed distribution of inhomogenities and anisotropics in the background 
radiation seems to be well described by the ACDM model, see e.g. [HE]- However 
they do not match perfectly with this model, [21 HI [S] , and consequently it is of interest 
to explore the complete phase-space of cosmological models and to see how differences 
from the standard picture affect the basic properties of the universe. 

Perturbations of anisotropic cosmological models have been considered by many 
authors, e.g., |H1 [Zl HI IB HH]: using methods depending on the choice of gauge, like 
the perturbation theory of Lifshitz and Khalatnikov [TT] or Bardeen's gauge invariant 
theory [12] . Unfortunately, the variables in Bardeen's theory are defined with respect 
to a particular coordinate system, making their geometrical and physical meaning 
unclear [13]. In the covariant approaches one circumvents these problems by using 
the spatial curvature rather than the metric as defining variables [Ml [15] and as 
the perturbed variables choosing objects that are zero on the background [HI fTZ] . 
and hence are gauge invariant [TB]. Some works on perturbations in anisotropic 
cosmological models along these lines can be found in [T9', . 

Bouncing cosmologies are of interest since the initial singularity can be avoided in 
these models. The observable universe seems to be close to homogeneous and isotropic, 
see e.g. |21| . Hence it may be approximated by the Friedmann-Lemaitre- Robertson- 
Walker (FLRW) metric that also serves as isotropic limits of different Bianchi models. 
In particular, the flat, open and closed FLRW universes are isotropic limits of Bianchi 
I, V and IX models, respectively |22]. From the Raychaudhuri equation it follows that 
a bounce is not possible in FLRW models if the strong energy condition is required 
to hold [23]. Nevertheless, when the weak energy condition holds a bounce can occur 
in a closed FLRW universe [24] , [25] , while matter violating the null energy condition 
can cause a bounce in flat or open universes 2B]-[1S]. Bouncing FLRW universes are 
also suggested by alternative theories of gravitation such as / (R) gravity [30] , / (T) 
gravity [31] and braneworlds '32], or by loop quantum gravity |33j-[36j. 

Cosmic microwave background (CMB) observations support a close to flat FLRW 
{\ftk\ < 10^^) universe [21]. Such constraints from the CMB observations have also 
been considered in bouncing cosmologies. Some models wherein a bounce appears due 
to a scalar field suffer from high tensor-to-scalar ratio [27], ^7\. However, this can be 
remedied both in general relativity ^38 and / (T) gravity [31| by introducing another 
massless scalar field. It was shown in loop quantum cosmology for a bounce followed 
by an infiationary phase that the tensor-to-scalar ratio places within the observational 
bound and that the low CMB multipoles are suppressed [35]- [H]. Actually, this 
suppression is observed but it can also have other origins Bouncing FLRW 

models with dust, radiation and cosmological constant are excluded by observations 
due to an elegant argument by Borner and Ehlers [68] . If the bounce took place before 
the formation of the quasars, that are observed at larger redshifts than 4, the present 
fraction of matter relative to the critical one would be less than 0.02, in contradiction 
with current estimates. 

It is not only the FLRW cosmologies that are interesting as possible models of the 
universe. The other models that have been investigated in some depth are those that 
have homogeneous but anisotropic geometries, namely the Bianchi and Kantowski- 
Sachs [43] cosmologies. Bianchi models were also used successfully to explain the 
CMB spectrum [11]-[17] and Kantowski-Sachs cosmologies might also be relevant for 
this purpose. The particular model presented in 4^ can explain some features of the 
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CMB spectrum, however not all. 

The cosmic no-hair conjecture says that spacetimes containing a positive 
cosmological constant finally evolves into de Sitter state. Wald proved that this 
happens for all Bianchi models with exception of Bianchi IX [33] ■ It was assumed 
in the proof that the fluid congruences are orthogonal to the homogeneous symmetry 
surfaces and that the matter energy-momentum tensor, not including the cosmological 
constant, satisfies the dominant and strong energy conditions. The isotropisation also 
happens under inflation provided via a scalar field with an exponential potential [50] - 
[55]. However, in the tilted case isotropisation does not necessarily take place even 
during inflation [24] , [53] . For Bianchi IX universes Wald gave a sufficient condition for 
them to evolve into de Sitter state. The evolution of Kantowski-Sachs universes with 
positive cosmological constant show similarities with that of Bianchi IX cosmologies 
[S3]. Although not the all initial conditions lead to de Sitter state, Moniz showed that 
the cosmic no-hair conjecture is widely valid in Kantowski-Sachs universes. Basically, 
the Kantowski-Sachs cosmologies with positive cosmological constant can evolve into 
de Sitter or Kasner states [24] . 

In this paper we investigate the density growth in perfect fluid Kantowski- 
Sachs cosmologies with positive cosmological constant mainly in those cases when 
isotropisation happens. We consider both bouncing and non-bouncing types of 
evolutions. The conditions under which bounces are possible were studied in [5S] 
and with a positive cosmological constant the Kantowski-Sachs models may under 
certain conditions undergo a bounce. This effect can also be achieved for Kantowski- 
Sachs solutions in i?" gravity, see [56]. To study density perturbations we use the 
1-1-3 and 1-1-1-1-2 covariant splits of spacetime [El [571 EHl US]- As inhomogeneity 
variables the spatial gradients of the density, the expansion, the shear scalar and one 
more auxiliary scalar to close the system, are used. These quantities are zero on the 
background, and are hence gauge invariant. By projecting along the preferred and 
orthogonal directions respectively, taking divergences of these projections and making 
harmonic decompositions of the spatial derivatives, the system is reduced to a first 
order system in time of four scalar quantities for each wave number. 

In section [3] we briefiy review the 1-1-3 and 1-1-1-1-2 covariant splits of spacetime. 
Then in section [3] the background solutions, of Kantowski-Sachs type, are studied 
and all background vacuum solutions are given. The perturbative equations are 
determined in section [4] and some analytical results are obtained in O Numerical 
studies of perturbations on different backgrounds, both with and without bounces, 
are performed in section jS] In [Appendix E| some results on bouncing closed FLRW 
universes are given for reference. 

2. The 1+3 and 1+1+2 covariant formalisms 

The 1+3 and 1+1+2 covariant splits of spacetime are here briefly reviewed. These 
formalisms are suitable also for perturbative calculations, as will be done in section|31 
For more details see, e.g., [16l[57] and [58l [59] , 

2.1. Preliminaries 

In |16[ 157] a covariant formalism for the 1+3 split of spacetimes with a preferred 
time- like vector, was developed. The projection operator onto the perpendicular 
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3-space is given by 

hi = gi + uau''. (1) 

Projections with of vectors are denoted by angle brackets = h^tp^ and the 

projected symmetric trace-free (PSTF) of a tensor is given by 



<ab> _ 



cd 



(2) 
(3) 



The covariant time derivative and projected derivative are given by 

V'a..6 = u''VcV'a...6 

and 

Vci^a...b = h{hi...htVf^d...e (4) 

respectively. The covariant derivative of the 4- velocity, w", can be decomposed as 



VqM;, = -Uailb + ^aUb = -Uailb + -^Okab + i^ab + Cfab 



(5) 



where iia = u^V^Ua is the acceleration, 9 = Vqw" the expansion, crab = V<aUb> the 
shear and LOab = [a^b] the vorticity of u"" . 

A formalism for a further split (1+2) with respect to a spatial vector n° (with 
u^-Ua = 0) was then developed in [58l[59]. Projections perpendicular to n° are made 
with 

Nl ^gl + Uau' - nan" . (6) 
Projected vectors can be decomposed with respect to n" as 

v<''> = Vn" + V (7) 

with 

V = n'''va and V = N''\i, = , (8) 

where a bar over an index denotes projection with N"-" . Similarly PSTF tensors ?A<af) 
can be decomposed as 



ib> 



where 



V'<a6> = * nam 



1 



^ab ^ 



':Na 



n°'n''il;^ab> 



2*(^n6) + ^ab , 



^<cd> ^ \^{ab} 



(9) 

(10) 
(11) 
(12) 



Here ^I*!"^! is symmetric and traceless. 

Derivatives along and perpendicular to n° are given by 

= n'VcV'a...6 = n''hlhi...hlV f^d...e (13) 

and 

Sc4'a...b = NiN^...N^Vf^Pd...e (14) 

respectively. Similarly to the derivative of the 4- velocity, the derivatives of Ua can be 
decomposed as 



Vanb = nattb + -^4>Nab + £,iab + Cab 



(15) 
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and 

ha = AUa + aa , (16) 

where 

aa = ha, 4'=6an°-, ^ = '^e°-''Sanb, (ab = S^ant], A — n^Ua, aa = ha (17) 
and 

Cafc = riabcn" = u'^rjdabcn" . (18) 
Here rjdabc is the totally anti-symmetric 4-dimensional volume element with 770123 — 

2.2. Fundamental equations in 1+3 split for the irrotational case 

The propagation and constraint equations are given in [57] . We will here only consider 
the case of a perfect fluid with vanishing vorticity. Imposing ujab = introduces only 
one new constraint 

ri-^'^'^VbUc = , (19) 

but with a barytropic equation of state, p = p(/i) where p is pressure and energy 
density in the rest frame of an observer, this equation is identically satisfied. From 
the Ricci identities one finds the following propagation equation for the expansion 

e - Vau" = -^02 + uau"- - 2(t2 - i(;t + 3p) + A , (20) 



where 



= ^a'^'aab (21) 

and A is the cosmological constant. 
The equation for the shear is 

^<ab> _ y<a^6> ^ ^Ig^jab i,<a^b> _ ^<ay>c _ ^ab ^ (22) 

3 

where Eab = Cacbdu'^u'^ is the electric part of the Weyl tensor. One also obtains the 
following constraints 

Vfecr"^ - ^V"6' = (23) 

H"'' = (curl aY^ = t^'"^<''\I , (24) 

where Hab = ^VadeC'^bcu'' is the magnetic part of the Weyl tensor. 
From the twice contracted Bianchi identities one obtains 

fi = -9iii+p) (25) 

VaP+(/i + p)ua = (26) 
and the remaining Bianchi identities give the following propagation equations 

£;«^''> = (curlH)"^ " + P)'^"'' ~ + 3cr<"c£^''>' + 2r]'"^<''uM^>d (27) 

ij<"''> + (curl Ef^ = -OH''^ + 3cr<%H''>'' - 2T^'"^<''ucE^>d (28) 
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and constraints 

S/,E-b - ivV - ii'^^'omH^'c = (29) 
Vfei/''^ + fi'^'^'aME^ = (30) 

respectively. 

By differentiating the constraints with respect to "time" and using the 
commutators between the "time" and "spatial" derivatives it was shown in [60] that 
the constraints are propagated for irrotational dust and in 61 this result was extended 
to the barytropic case p — p{n)- 



3. Kantowski-Sachs 

Kantowski-Sachs cosmologies, [13], have a 4-dimensional isometry group acting 
multiply transitive on 3-spaces with topology i? x 5*2, i.e. they are locally rotationally 
symmetric (LRS). With zero vorticity the line-element can be written as 

ds^ = -dt^ + al{t)dz^ + al{t) {dd^ + sir? M(p^) (31) 

with the 4-velocity of matter given by w = ^ and the direction of anisotropy by 
n = ^ ^ . The expansion is given by 

+ (32) 
ai 02 

and in the tetrad 

uj^ — dt, oj^ ~ aidz, lJ^ — aidd^ — ai sini^di^, (33) 
the shear is given by 

E ^ a„ = -2^22 = -2a33 = ^ f ^ - ^ V (34) 
3.1. The evolution equations 

Due to the LRS symmetry the shear and electric part of the Weyl tensor can be written 
as 

aab = ^{UaUt - ^Nab) ■ (35) 

and 

Eab = S {UaUb - ^Nab) (36) 

respectively in terms of the anisotropy vector and the projection operator Nab- 
Given an equation of state and the cosmological constant, the Kantowski-Sachs models 
are completely determined in terms of shear, E, expansion, 9, and energy density, /i. 
The electric part of the Weyl tensor is given algebraically as 

£= -^ti-^A-^' + ^e' + l^e, (37) 

whereas 

Hab ^ Ua = aa = (f) = (. = Cab = A = aa = . (38) 
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From the equations in section 12.21 the fohowing evohition equations are obtained 

t = p, + h + ^j:^-j:e-^6' (39) 

fl=-9{^I+p) (40) 
Alternatively one of the equations can be replaced by 

k 

where K, given by 



l(^_^3p_2A)_^I]2. (41) 



k^-{'^9-Y)K (42) 



K = h + A+^J:^ -^9'^ >0, (43) 



is the curvature of the 2-spheres 82- 



3.2. Dynamical system analysis of background 



Dynamical system analysis of the Kantowski-Sachs models with positive cosmological 
constant were done in, e.g. [62j[24]. We here follow the notation in [24]. The relations 
between their variables and ours are given by 



^ - ^ o - ^ 
The equilibrium points are given by 



Q+= -T^^ ^^-^^ = (44) 



±F: Flat Friedmann: A = S = 0, /i = 6*^/3, Qo — ±1, saddle points 
+K±: Kasner: A = = 0, S = t|^, Qo — 1; Q+ — ilj source points 
-K±: Kasner: A = /i = 0, E = ^|^' Qo — ~1: Q+ = ±li sink points 
±^5: de Sitter: 11 — T, ^ 0, 9 ~ ±\/3A, sink/source points 
±X: fi = 0, 9 = ±\/A, E = ±|\/A, saddle points. 

The points ±X are exact vacuum solutions with metrics 1 241 



ds^ = -dt^ + e^^^^dz^ + ^ (rfi?2 + sin^ Mip^) . (45) 



3.3. Vacuum solutions 



In the vacuum case the equations can be integrated completely. If the curvature of 
S2, K, is a constant it follows that E = (or X = 0). The system (I3n])-(|31]) then 
reduces to 9 = A — 9^ with solutions 

9 = VA^^^ , (46) 



where C is a constant of integration (and 9 — \/ A, corresponding to the critical point 
-i-X). By redefining the time coordinate t ^ t + to, these solutions can be rewritten 
as one of the following two 



cosh(\/At) sinh(\/Ai) 
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(C = gives the critical point -X with 9 — — \/A.) The corresponding hne-element is 
given by 

ds^ = -dt^ + f{t)dz^ + ^ {dd^ + sin^ Mtp^) . (48) 

where f{t) for the first sohition is given by 

f{t) = ao cosh(%/At) . (49) 

These spacetimes experience a bounce in the z-dirction at i = (and are non- 
expanding in the perpendicular directions). They start at the critical point _X and 
end at+X. For the other solution one gets 

f{t) = ao sinh(\/At) (50) 

and hence these solutions are singular for t — Q. They start at the critical point 
(Kasner) and ends at or (for large negative t) start at _X and end at 

If K is not a constant it can be used as the independent coordinate. The system 
(|39 p - p2|) for /i = p = 0, together with the constraint p3|) . then reduces to (see [63] 
for details) 

^ _ _^ 2A - 4i^ + QMK^I'^ S - ± ^^^^'^ f51) 

2^2MK^'^ - K + k/?,' ~ Z^2MK^/^ - K + K/Z ' 

where M is a constant of integration. Changing independent coordinate through 
K = the corresponding line-element becomes 

If A = 2M/T — 1 + > this is the space-homogeneous region of the 

Schwarzschild-de Sitter metric [B^, and M = gives de Sitter space. With 
AM^/3 > 1/27 the requirement A > is satisfied for all positive T if M > and the 
metric starts at the critical point +1^+ (T — 0) and ends at -^dS (de Sitter) . A > 
also in the region {~oo,Tq), where Tq < is the only real solution of A = 0. These 
spacetimes start at ^dS and ends at If Af < the situation is reversed so that 

for T < there is a class of solutions starting at _ dS for large negative T and ending 
at -K- for T = 0, and another class starting at +A'_ for a positive Tq and ending at 
+^5* for large positive T. 

For < AAf^/3 < 1/27 and M > there is one class of solutions starting at 
-dS at negative infinity and ending at -Kj^. for a negative Tq, another class starting 
at +K+ for T = and ending at -K+ for a positive Ti and a third class starting at 
for a positive T2 and ending at -^dS for large positive T. Similarly, for negative 
M, there is one class going imm^dS to^K+, a second from+_ft'_ to^K^ and a third 
from+ii'_ to ^dS. 



3.4- Solutions with matter 

The general dust solutions in terms of elliptic functions were found in [BSl [SB]. A 
few of them can be given in terms of elementary functions, see e.g. [661 167] . For 
certain choices of the parameters in [66] physically reasonably solutions starting from 
a pancake like singularity can be obtained. These solutions end in expanding de Sitter. 
The solution [67] starts in +F and ends in -^X. 
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We here consider approximate solutions with /i ^ p ^ A. To first order the 
dependent quantities are written as 9 = 9q + 6i, E = Eg + ^^i a-nd /i = /Zi. Assuming 
a hnear equation of state p = (7 — and Eq — |0o for the background vacuum 
metric, the foUowing first order system 

Ai = - 0o7Mi (53) 

El - ^^1 = 7/11 + 00 (^Si - ^01^ (54) 

Si + ^01 = i(6 - 57)mi - 200 (^Si + ^01^ (55) 

is obtained. In appendix [Appendix A| solutions around the bounce metric (j49p and 
the critical points ± AT, (|^5|) . are given. Since for all of these solutions some terms grow 
unbounded, they are only valid for limited time intervals. Nevertheless they can be 
used to check our numerical codes for shorter time intervals in the case with densities 
that are initially low, and also for finding suitable starting conditions. 

3.5. Bouncing solutions 

Bouncing FLRW universes containing dust and radiation with cosmological constant 
are excluded by observations [55] , see also [Appendix E[ This relies on a very simple 
argument. If quasars are observed at redshift z > 4, the dust component fl„i cannot 
overcome the value 0.02, otherwise the bounce would occur at a smaller redshift than 4. 
We now consider whether a similar argument also holds for bouncing Kantowski-Sachs 
universes with cosmological constant. The bounce occurs either 

i) in the direction of anisotropy at ai* for which di* = and di* > 
or 

ii) in the perpendicular direction at 02* for which 0,2* = and 02* > 0. 

By defining the average scale factor a through 

= 3-. (56) 
a 

and using p2l) one obtains — aia^. If we assume that the universe is filled with 
dust (7 = 1) and radiation (7 ~ 4/3) Eq. is satisfied by the following density 

Ai = Mm + /^r = MmO (j^^ + MrO (j^^ , (57) 

where the subscript denotes present values. We then introduce the usual 
dimensionless parameters 

f^m = 3/i,„o/0o , = 3Mro/0o and ilA = 3A/02 . (58) 

Case i): Since 0/3 + E = di/ai, we have < 0*/3 + E* and by using Eqs. dSHl), dH]) 
and the notations (l58l) this leads to 



Eq. (US]) gives 



3 di* r2,„ / flo \ ^ r^r- / oo 

00 ai* 2 \a* / 3 \a. 



0<^— = ^(-) +^(-) +^^A- (59) 
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and another combination of Eqs. (P^l) and (HTj) . 29^/3 — E*, gives 

-^'^ = nj'-^)\inf!^)\ (61) 

The inequahties are always satisfied with a positive cosmological constant and the 
equations also give the values of 0,2*, 02* and ai* at the bounce. Therefore, we have 
no other constraints unless we can find some further integrals to the system ([55)1 - ([1T|) . 

We can also consider the equations (|39|) . (|4T|) and p3)) at the present time to- 
Neglecting radiation one obtains 



-qiO = -ilm + iiA (62) 



2H10H20 



-3920 = g^2 (63) 



1 i IS.H 1 

— 5-^=a„ + f7A-l + ^-^, (64) 

where i/o = clq/oq is the present value of the Hubble constant, Hiq = ctio/aio, 
H20 = (120/0,20, AHq = Hio - H20, qio = -aio/{aioH^) and ^20 = -020/(020^^0)- 
How the quantities Hio, H20, Qio and 520 relate to observations can be seen from the 
general expression for redshift (when both emitter (E) and receiver (i?) have fixed 
spatial coordinates) 



l + (65) 

(see for example [W|)- Here Pa = gabdx^/dX (with affine parameter A) denotes the 
covariant components of photon 4-momentum. By integrating the geodesic equations 
for the Kantowski-Sachs metric the following redshift formula is obtained for photons 
moving in the 9 = t:/2 plane 

„/o r 9 . 9 91^/2 
ao faoy''^ 



l + z^ — , — I 

a \ a / 



(66) 



where a = ai/a2 and Pz and p^p are integration constants. For example, photons 
moving in the z-direction have a redshift 1 + z = oio/oi whereas those moving along 
the 2-spheres have redshifts 1 + z = 020 / 0-2 ■ 

In the isotropic limit AiJo 0, A^/o = <Zio — 920 — ^ and l/a^o — the 
equations ((62)) - (|64)) agree with those of the flat FLRW model with cosmological 
constant (go = 0000/^0 ~ ^rn/2 — JIa, 1 = ^rn + ^^a)- Since there are indications 
of anistropies, but with large uncertainties, both in the Hubble and deceleration 
parameters, see e.g. [TOj and [7T], it would be of interest to get better estimates 
of these parameters, even if present values do not seem to support Kantowski-Sachs 
models. 

Case a): From K — l/a\ we find K — 2 (—02/02 + ia^/a^ K and therefore Ki, < 0. 
Then, taking the derivative of Eq. ([^^ to get another expression for K and using 
Eqs. (IMl) and (gT]) together with 9^, = 3E*/2, we find 

nm<0 . (67) 

Hence a bounce in the scale factor 02 is not possible, whereas a bounce in the direction 
of anistropy cannot be excluded simply from this type of argument. 



Density growth in Kantowski-Sachs cosmologies with cosmological constant 11 

If we instead look at a bounce in the average scale factor at a = a^, so that 
9^ = ia^/a^ — and > 0, the following inequality 

2 1 r. 3^ /ao^^ 



can be derived, but no upper bound for the redshift at the bounce can be obtained 
from this either. 

4. Perturbations on Kantowski-Sachs 

In this section we calculate the equations governing the growth of density perturbations 
on a Kantowski-Sachs background to first order. The inhomogeneities will be described 
by quantities that are zero on the background, and hence are gauge invariant [18) . The 
primary variable is the density gradient 

Pa ^ ^ . (69) 

M 

Here a is the average scale factor, defined in (1551) . The density fluctuations ^ 
on a length scale I are related to the quantity Va through ^ ^ [VaV^Y^^l / a — 
{VaVfl'^lo, where Iq = l/a is the comoving dimensionless length scale. However, 
note that the quantity 

Si,) ^ (70) 

depends on the identification between the fictitious background with density and 
the real universe with density and can be given given any value by changing the 
identification [57], whereas Va is gauge invariant. 

To close the system three auxiliary quantities, that we choose as 

Za = aVad, Ta = aVa<J^ and 5a = aVaK'^af,) (71) 
will be needed. Here the traceless part of the 3-Ricci tensor is given by 

Sab = '^Rab — 2 ^Rhab = —(^<ab> — 0<^ab + ^ <aUb> + U<aUb> 

= Eab + <y<a<yb>c — -^daab ■ (72) 

In accordance with ([9|) aab and Sab can be decomposed as 

aab = ^{riaUb - ^Nab) + '^■^(anb) + Sab (73) 

and 

Sab = SiUaUb - ^Nab) + "^-S^anb) + Sab (74) 

respectively. Note that Sab to zeroth order is given in terms of the other quantities, 
but to first order is an independent quantity. 
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4.1. The first order equations 

The propagation equations for the gradients are obtained by taking the gradients Vq 
of the propagation equations in section 12.21 and then using the commutator between 
"time" and "spatial" derivatives acting on a scalar, that to first order reduces to [57] 

Va(/) - (Va/)- = -Uaf + ^O^af + <Ja cf ■ (75) 

In [Appendix B| some details of the calculations are given. Finally the following first 
order system 

Va = —Va - l^UaU'V, + ^SPa - 2,(1 + ^) (76) 

/i 2 2 /i 

2z,-T l..-T^ , 1^^ , 3 pip' f~ ^ 



3 2 z 2 2/i + p\2/ 

-VaV''2?b (77) 



M + P 



Ta = -2eTa - l^naU-Tc + ^ETa ~ Sa ~ ^J^^ Za + ( S + -^2?„ 

2 2 2 2\ / fi + p 



S(-n,n''V,V^Pfc - -VaV'Pb) (78) 



'2 " 2 

5„= (E2+2|Jjra+ (^E-^0-2|j5a-^Enan^5,-s('^5 + ^E0)Z, 



lj,+p \2 2 2 y " ^ " n + p 



where 



^En''n=V„VbZe - ^EVaV^Zb - V^V'-Tb , (79) 



^ = --^--A-iE2 + -02^--i^<O (80) 
3^ 3 2 9 3 ^ ^ 

to zeroth order and p' = dp/dfi, is obtained. As seen two apparently singular terms, 

28"^ /Y?Ta and — 2S'/E5a, appear in ([79)) . Near points where E = it is hence suitable 

to remove these terms by changing the dependent variable 7^ to 7^ through 

Ta - E^ta + . (81) 

Instead factors 1/5 and 1/5^ will be introduced, but since S = —2K/'i < these will 
be well behaved. 

4-2. The projected equations 

The equations ([76| - (|79l) can be decomposed into two sets by projecting with Ua and 
Nab respectively. For Kantowski-Sachs, for which (|38l) holds, it follows from equations 
(|16p and ([Tsl) for the derivatives ha and Va".6 of Ua that 

ha = VaUb = . (82) 

Since = and habu"" = it also follows that the derivatives of A^ab, Nab and WcNab, 
become zero: 

A^ab = U^Vcigab + UaUb + ria^b) = UaUb + U&Wa = 
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VcA^afc = Vc(mqU6) = hJh^h^{udV fu^ + u^V fU^) = 0. 

Hence, since we only need Ua and Nab to zeroth order, we can just let and Nab 
"pass through" the derivatives when projecting equations ([75)1 - ([75)1 . 
With the definitions 

V = Van'' , Z = Zan" , T = Tan"" , S = San"- (83) 

the equations projected along n° become: 



Z= - \^V-(le + Y\z-2T+l^(~S+lYAv-J^n-VaV'Vb 



2 V3 / + 2 / ^i+p 



r = - (26* + S) r - -S"Z + --^S + 6IS I? - 5 - -^E(-7i''Van''V6l? - 
2 2/i + p\ / /i+p2 



£2/ I 3 E/ V2 Sy ^i + p \2 2 2 



The terms V^Vb etc. can be decomposed as 
to first order. 

The orthogonal equations, obtained by projecting with Nab, can be found in 
[Appendix C[ 

4-3. Scalar equations 

To treat the spatial derivatives appearing in the equations we will do a harmonic 
decomposition. For this purpose it is suitable to get the spatial derivatives in the 
form of two Laplace-like operators, 

^2 = SaS" and A = n^VaU^^b (89) 

acting on our variables. To obtain this we define new variables 

V = n^VaV and p = S^Va (90) 

and similarly for the other variables. We then act on the system ([M)) -(|57 | with the 
operator n^Va and use the commutation relation 

*-*=Q0 + S^*, (91) 

[5U] . that holds to first order. To remove the singular terms 2S'^ /T?T and — 2S'/I]5 
we now also make the aforementioned change of the dependent variables T and y to 
T and "Z" through 

f=T?f+^S and r = Y?f+^S (92) 
S S 
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respectively. The system for the hat variables then becomes 



V - 

Z : 
f 





{-- 








e 




- 25] 















1 



^3 
MP' 1 



21]) Z - 

+ 2EH 
1. 



1 



7 A* 



3 ^lp' 



IJ- + PS 
^A(5 + ^) 



S3 
S 

-5]-i 
2 3 



r- 

] AV 



5 + iE> 



P- 2^5- 2S2f 



5 



A(^) 



fJ-P 



M + P 



■A 



1 



v + p 



3T 



|A(t + f) 



(93) 
(94) 

i + 

(95) 



MP 
M+P' 



5( :loi: + Is-Iy^ 



25*2 r 



S - 26*") 5 + EA ("i - ) - I]2 A ( r + T 



5 
1 



MP 
M + P 



-5 + ieis - -Y? 

3 2 



AX> 



1 



- - -6*1] + 2E2 A (^) 



|a(5 + ^). 



(96) 



By taking the 2-divergence of the system (|C.ip - (jC.4[) and using the the commutation 
relation 

'1 „ 1. 



(97) 



[5U] . (where (5°5'a = ^ according to the above definition) a similar system, that can 
be found in [Appendix C[ is obtained for the slashed variables. As we will see in the 
next section, the hat and slash variables are closely related. 

Since the scale factors ai{t) and 02 (i) appear in the spatial derivatives, the time 
dependence of the variables V and p will go as V/ai and I?a/a2 respectively. This is 
most easily seen by calculating P and p in the tetrad (|33l) . giving 

1 dD-^' 



V= — and p = — ' 2 

ai oz 02 



dd sin I? dif 



Hence the time dependences of the variables 
2?j| = aiD and — a2p 



(98) 
(99) 



give a better representation of the development of the relative density contrast ^ . 



4-. 4-. Harmonic decomposition 

We will use a harmonic decomposition 



fcii ,k± 



of the dependent variables. Here satisfies 

fc2 

APfc„=-4Pfc„, 5aPk,^Pk,=Q 



(100) 



(101) 



Density growth in Kantowski-Sachs cosmologies with cosmological constant 



15 



where fcy are the constant co-moving wave numbers in the direction of anisotropy and 
the scale factor in this direction, ai, is given by the zeroth order equation 

— ^iez + E. (102) 
ai 3 

If a coordinate z is adopted to the 1-direction a possible choice for Pk^^ is Pk^^ = e**^"^ 
as can be seen by direct substitution into (jlOll) . using a tetrad, e.g. ((33)) . adopted to 
the symmetries. Similarly the harmonics Qk^ are introduced on the 2-sheets as was 
done in ^P]. They satisfy 

1.2 

6^Qk^ = — iQfc^ , Qk^ = Qk^ = (103) 
where the scale factor on the 2-sheets, 02, is obtained from 

^ = ^^-^S, (104) 
02 3 2 

and fcj_ are the co-moving wave numbers in the perpendicular directions. 

In the irrotational case there is a simple relation between the coefficients of the 
modes for the hat and slash scalars. We note that to first order V and p can be 
written as 

V = -An and p = -S^u (105) 

respectively, with similar expressions for the other scalars. From the commutation 
relations in [59] it follows that the operators A and commute to first order when 
the vorticity is zero. Hence 

S^t) - 6^ (^^A^i^ = A - (106) 

holds to first order. Using the harmonic decomposition (|100|) and equations (jlOip and 
(fT03l) . one has 

,2 fc2 
Go Oi 

with similar expressions for the other scalars. Hence we find the following relations 

2 / X 2 



where ^ is P, Z, T or 5 and ^ is ^, ^, /" or between the coefficients of the 
different modes. In the following we will often suppress the subscripts fc|| k±^ when it 
is obvious that we are refering to harmonic coefficients. 



5. Perturbations of analytical solutions 



In this section we summarize some results on perturbations around a few of the 
exact vacuum solutions from section 13.31 in the infinite wavelength {k — 0) limit. 
More details are given in appendix [Appendix D[ It might seem unphysical to consider 
density perturbations in a vacuum solution, but the solutions might approximate 
perturbations around some background metrics with p <C /i <C A. Furthermore, the 
analytical results are useful for a comparision with results of the numerical codes. This 
also gives a method of identifying suitable initial conditions. 
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5.0.1. Perturbations of With the +X-solution as background one finds the 
fohowing fourth order equation for the density gradients 2? in the direction of 
anisotropy 

+ ^ VAP(3) + ^Ai> + — A3/2p + ^A^T> = (109) 
3 3 27 81 

if the harmonic numbers are put to zero (corresponding to large wavelenghts of the 

perturbations). The solution , i.e. V muhiphed with the scale factor ai cx e^*, is 

given by 

- Aie-8v^*/3 _^ ^ ^^^^ ^ A4)e-2v^*/3^ (110) 

where Ai, A2, A3 and ^4 are integration constants (or slowly varying functions over 
space if the wave number is not exactly zero). Hence there is one growing and three 
decaying modes. 

The solution for the density gradients 2?j^ = in the perpendicular directions 
is then obtained from V±_ = as V cx e^^X* as 

Vi_ = Sie-5^*/3 + B2e4^*/3 ^ (^^^ ^ Bi)e^'/\ (111) 
As seen, in these directions that are not expanding, the modes are growing faster. 

5.0.2. Perturbations of vacuum bounce solution With the vacuum bounce solution 
(|49p as background it is suitable to change the independent variable to 6 through 

The following equation for the density gradient V in the direction of anisotropy 
d^V (76(2 + 6A) d^p (704 - 6A6l2 - 6A2) 5{8e^ - ibAO'^ + 18A'^) dV 

'W ^ 36i(e'2 - A) rf^ 36»2(6»2 _ A)2 27e{e'^ - A)3 

(-735A6'4 + 495A26'2 - 406*6 - 270A3), 



8102(5»2 _ A) 



4 



-X> = (112) 



is then obtained in the long wavelength limit. The solution — aiD, where 
ai cx l/VA — is then given by 

- (Ai + ^20) (A - 02)1/3 ^ _ ^2) -1/6 + (A - 02)1/3 X 



1 / 02 

In 1 




2 \ A) 4VA V \/A - / VA - 02 VVA 



(113) 



(assuming 02 < A). Here = \/A tanh(\/At) and hence A - 02 = A/ cosh2(\/At). 
The mode Ai starts growing, obtains its largest value at the bounce and then starts 
decaying. The A2 mode also starts growing, but reaches its maximum before the 
bounce. It then decays to zero magnitude at the bounce and after this passes through 
a new maximum before it eventually decays. The A^ and Ai modes initially decay, 
pass through zero at i = and then grow unboundedly. 

As in the previous case the growth of the density perturbations in the non- 
expanding directions are obtained from 2?j_ = a2p = 02 (f^) ^/(A — 0^) 
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as 

2?± = (Si + B29) (A - -1/6 + BsO (A - 6^^^^ + B4 (A - 

- In 1 - — ^ In ^ + , arcsm —= . (114) 

2 V a; 4Va vVa-^^; va^ i m v y 

These modes decay before the bounce and grow after. The modes B3 and B4 are close 
to zero for a fonger perfod of time around the bounce. 

For this solution perturbations with co- moving wave number fc|| > 2ai(io)VA 
cross the horizon in the anisotropy direction within a finite time, and one would 
expect a difference in behaviour between perturbations with wave numbers that are 
smaller or larger than 2ai(to)\/A, respectively. This is also favored by the numerical 
analysis in section [51 where typically the perturbations in similar models show an 
oscillatory behaviour for wave numbers larger than this value. 



6. Numerical solutions 



The time evolution of the density gradients are solved for numerically in some 
representative cases. To find suitable backgrounds we have used the phase space 
analysis in [5J (see also section IX^ . Source points are given by the expanding Kasner 
metrics +K± and contracting de Sitter -^dS, and sink points by the contracting Kasner 
metrics -K± and expanding de Sitter -^dS. The plots shown below either start from 
expanding Kasner +K±, or orginates from them. They mainly end at expanding de 
Sitter +dS, being a more likely state of the late universe, but we also give examples 
of metrics starting from +K- and ending in -K-. Saddle points are given by the fiat 
Friedmann models ±F and the solutions ±X, which were useful in finding solutions 
undergoing a bounce. This is related to that the vacuum bounce solution (|49p starts 
in _X and ends at ^X. For each type of path in phase space one case with radiation 
and one with dust are given. 

The system contains four different modes, Z, T and S, for each background 
and wave number. In most of the presented cases just an initial density perturbation 
is assumed, but in the two first cases we also show the effects of an initial pure 
shear perturbation. The behaviour for different types of initial perturbations has 
been checked, but the differences one can see are of the same types as those given 
between figures [2] and [H and figures |4] and [5] respectively. The differences are largest 
in the dust case, where the density perturbations depend more strongly on wave 
number with initial shear and expansion perturbations than with an intial density 

perturbation due to that factors like A + 'JT^ appear as source terms in the 
equations. For each radiation case we show the evolution for the wave numbers 
k = A:||/aio = k±/ci20 = 0,1,5 and 20, including both sub horizontal and super 
horizontal perturbations. For dust, where the dependence on wave number is small 
for intial pure density perturbations, we only show the evoultion for one wave number 
except for the case with intial shear perturbations. 

Throughout this section we choose the cosmological constant A = 1 . Other initial 
values are given in the figure captions or text. The code was tested against the results 
of section [5] 



6.1. Solutions with a bounce 
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Figure 1. A radiation background that experiences a bounce. Initially it is close 
to _X and a late times it approaches +dS. Initial values at to = 1 are given by 
HO = 0.02, 6*0 = -0.6, So = -0.6, e» = -0.8 and = 0.1. 




Figure 2. The growth of the density perturbations in the background given 
by figure[l]for the wave numbers k = fcy/aio = kj_/a2o = 0, 1, 5 and 20. Initially, 
at to = I, {> = 'p = 0.001. The time of the bounce is indicated with a dotted 
vertical line. (Here aio = ai(to) and a20 = a2(to).) 
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Figure 3. The growth of the density perturbations in the background given 
by figure[l]for the wave numbers k = fcy/aio = fcx/i20 = 0, 1, 5 and 20. Initially, 

at to = 1, T = X = 0.001. The time of the bounce is indicated with a dotted 
vertical line. 

6.1.1. From _X via j^X to +dS, radiation This background, where the equation of 
state is given by p = /i/3 (radiation), passes through a state where the directions 
perpendicular to the direction of anisotropy have small and decreasing expansion 
rates, 9± = 02/02, that become almost neglible for a period of time. After this 
the expansion starts again to finally reach a constant value 9±_ — -y/A/S. The 
anisotropy direction comes from a contracting state, goes through a bounce and then 
starts expanding. Eventually the expansion rate in this direction, 6\\ = di/ai, also 
approaches 6'|[ = y^A/3, giving a total expansion rate of = VSA. The initial state 
is close to the critical point -X and for an intermediate period the solution is close 
to +A, but for large times the solution approaches the sink point de Sitter -^dS. The 
evolutions of density, fj,o, expansion, 9q, expansion in the anisotropy direction, ^y, 
expansion in one of the perpendicular directions, 9j_ and shear Eq are depicted in 
figure [TJ Note that, as is found by running time backwards in the numerical code, 
that this solution can be seen as originating from the source point +K^. 

In figure[2]the growth of density perturbations for different values of the co- moving 

wave numbers fcy and k± is shown. Since from equation (|108p V± = (^^^ a^'^W 
only show . The initial values of the density perturbations at to = 1 are given by 

V = p ^ 0.001 and the other quantites are initially put to zero, Z = = T ^ 
S = 1^ = 0. For A; = the density gradient in the direction of anisotropy reaches 
a small maximum after the bounce and after this a small minimum before it starts 
growing unboundedly. In the perpendicular directions the gradient for the fc = case 
is roughly constant before it also starts growing. For higher values of the wave number 
k the density gradient in the anisotropy direction shows an oscillatory behaviour with 
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Figure 4. The growth of the density perturbations "Dn and in the dust 
background of section 16.1.21 for the wave numbers k = fcy/aio = fcx/120 = 1- 
Initially, a.t to = 1, V = p = 0.001 and ^to = 0.02, 6*0 = -0.6, Sq = -0.6, 
6|| = —0.8 and = 0.1. The time of the bounce is indicated with a dotted 
vertical line. 

an initally increasing amplitude that later on decreases, but does not fall of to zero. 
This behaviour with a local maximum in the density gradient at or slightly after the 
bounce in the bouncing direction seems to be typical. For corresponding k values in 
the perpendicular directions the oscillations initially have an approximately constant 
amplitude that with time slowly starts growing. 

In [Appendix E| the corresponding growth of density perturbations in a closed and 
isotropic radiation filled universe undergoing a bounce (now in all three directions) is 
shown in figure [E2l The local extrema seen in figure [2] are absent in the isotropic case, 
apart from a small dip for the fc = mode, whereas the behaviours for large times, 
when the universes in both cases approach de Sitter spacetime, are similar. 

In figure [3] the growth of density perturbations for a case where the density 
perturbations initially are zero, = ^ = 0, is shown. Instead there is an initial 

shear perturbation T = 7" — 0.001 (and Z = ^ = S^^ = 0). As seen the super 
horizon modes fc = 0, 1 grows unboundedly, whereas for higher fc- values the amplitudes 
remain small. 

6.1.2. From via j^X to +(13, dust The conditions here are similar to those in 
the previous section, but the equation of state is now given by dust, p = 0. With the 
same initial conditions for the background quantities as in the radiation case, their 
evolutions are close to those of the radiation case, see figure [TJ First we show the time 
development of the density perturbations for the case when initially T) = 'p = 0.001 

and Z = .^ = T = /' = 5 = <^ = 0. Since in this case the evolution of the density 
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Figure 5. The growth of the density perturbations in the dust background 
of section 16.1.21 for the wave numbers k = fc||/aio = fc_L/a20 = 0,1,5 and 20. 

Initially, at to = 1, 7^ = 7" = 0.001. The time of the bounce is indicated with a 
dotted vertical line. 

perturbations is relatively insensitive to wave number, we here only give them for k = 1 
in figure 31 In the direction of anisotropy the density gradients first grow and reaches 
a maximum at approximately the time of the bounce and then decay into a small 
but nonzero value. For the perpendicular directions the behaviour is inversed. First 
a minimum is reached around the time of the bounce, and then the gradient grows 
towards a constant value. A comparison with the evolution of density perturbations 
in a bouncing closed and isotropic dust universe, see figure lE3l in [Appendix E} shows 
that the local minimum (maximum) is absent in the isotropic case. 

Next, in picture [SI the growth of the density perturbations for the case when 

there is an initial shear perturbation T = '~f = 0.001 (and the other perturbations are 
zero), is shown. In this case the density perturbations depend on the wave number 
k. The reason for this can be seen from the systems of equations (|93l) - (|96l) and (|C.5|) - 

(jC.Sp . where A {t + 7") and 6"^ (t + /"^ appear as source terms. 

Note that all perturbations asymptotically approach constant values. This is 
consistent with that the background asymptotically approaches de Sitter. The result 
holds also for the cases 16.3.21 and I6.3.4[ that both are dust solutions approaching de 
Sitter. That this result does not apply for the radiation solutions can be understood 
from that terms like p/ fJ- = 1/3 remain in the equations also in the limit /i — > 0. 

6.2. Recollapsing solutions 

6.2.1. From to -K^, radiation Here the background, that is filled with 

radiation, starts close to an expanding Kasner, +/^-, reaches a largest value, and 
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Figure 6. A radiation background that starts expanding, reaciies a largest value 
and then recoUapses. It starts close to +K- and ends at -K—. Initial values at 
to = 1 Eife given by fio = 0.2, 9o = 3, So = 2, 8a = 3 and 9± = 0. 
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Figure 7. The growth of the density perturbations in the background given 
by figure[6]for the wave numbers k = fcy/aio = k±/a2Q = 0, 1, 5 and 20. Initially 
ti = p = 0.001. 
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Figure 8. The growth of density perturbations 'D^^ and in th^ dust 

background of section 16.2.21 for the wave numbers k = fcy/aio = fcx/120 = 1- 
Initially V = p = 0.001 and /xq = 0.2, 60 = 3, Sq = 2, ey = 3 and 6»x = 0. 

then approaches a coUapsing Kasner, In the direction of anisotropy it is always 

expanding, whereas it is collapsing in the perpendicular directions. Initially this 
collapse velocity is negligible, but eventually it will dominate over the expansion in 
the direction of anisotropy. The evolutions of density, expansion, expansion 
in anisotropy direction, , expansion in one of the perpendicular directions, 9± and 
shear Eq are depicted in figure [51 In figure [7] the growth of the density perturbations 
for different values of the comoving wave numbers /c|| and k± is shown. Initially 

V = 0.001 andZ = ^ = f = 7' = 5 = <^ = 0. 

As can be seen the density gradient in the direction of anisotropy decays for all 
fc-numbers. In the collapsing directions the gradients grow. 

6.2.2. From +A'_ to -K^, dust This is a similar situation to the previous one, but 
with dust. The evolutions of the background quantities are once again close to those 
of the corresponding radiation case in 16.2.11 see figure [51 Also in this dust case the 
perturbations are rather insensitive to the wave number, and we hence only give them 
for one /c-value in figure [51 As before the the density perturbations decay in the 
expanding direction and initially grow in the contracting directions. 

6.3. Solutions without a bounce 

For all cases in this section the initial conditions of the perturbations are "D = Tp = 
0.001 and i = ^ = r = T = 5 = <^ = 0. 
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Figure 10. The growth of the density perturbations in the background given 
by figure[9]for the wave numbers k = fcy/aio = k±/a2Q = 0, 1, 5 and 20. Initially 
■D = -p = 0.001. 
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Figure 11. The growth of density perturbations TJ^^ and in the dust 

background of section 16.3.21 for the wave numbers k = fcy/aio = fcx/120 = 1- 
Initially V = p = 0.001 and /^o = 0.36, 6»o = 0.6 and Eq = 0. 

6.3.1. Friedmann to de Sitter, radiation This radiation background passes close to 
an expanding Friedmann, -^F, and ends at an expanding de Sitter, -j-dS. It originates 
from an expanding Kasner, as can be seen by integrating backwards in time 
with the same initial coditions. Note also that at t = 1, the starting point in figure [HI 
the growth rate of the shear is nonzero. The evolutions of density, /iq, expansion, Oq, 
expansion in anisotropy direction, , expansion in one of the perpendicular directions, 
9± and shear Sq are depicted in figure |9l The expansion in the direction of anisotropy 
is for a period of time dominating over the expansion in the perpendicular directions. 
In figure [TO] the growth of the density perturbations for different values of the 
comoving wave numbers fcy and k± is shown. In the direction of anistropy the super 
horizon modes fc = 0, 1 eventually grows unboundedly, whereas for higher modes the 
amplitude slowly falls of. In the perpendicular directions the fc = 0, 1 modes grows 
faster and the amplitudes of the higher modes slowly grow. 

6.3.2. Friedmann to de Sitter, dust For the corresponding dust background to the 
radiation case of section 16.3.11 the evolutions of the background quantities are once 
again close to those of the radiation case, see figure [9] As in the previous dust models 
the perturbations are insensitive to the wave number. In figure [TT] the growth of 
density perturbations for the comoving wave number k = 1 are shown. The density 
gradients in the direction of anisotropy falls of to a nonzero constant value, whereas 
in the perpendicular directions they grow towards a higher constant value. 

6.3.3. +K- to +dS, radiation This radiation background starts close to an 
expanding Kasner, +K^, and ends at an expanding de Sitter, -^dS. The evolutions of 
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Figure 12. A radiation background starting close to expanding Kasner, -\-K—, 
and ending at expanding de Sitter +dS. Initial values at to = 1 are given by 
^10 = 0.2, 6*0 = 3.1 and S = 2. 
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Figure 13. The growth of the density perturbations in the background given 
by figure [T2I for the wave numbers k = fc||/aio = fcx/ci20 = 0, 1, 5 and 20. Initially 
■D = p = 0.001. 
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Figure 14. The growth of density perturbations X>|| and X>x in the dust 
background in section 16.3.41 for the wave numbers k = A;||/aio = k±/a20 = 1. 
Initially V = p = 0.001 and /xq = 0.2, = 3.1 and E = 2. 

density, /xq, expansion, 9o, expansion in anisotropy direction, expansion in one of 
the perpendicular directions, 9± and shear Sq are depicted in figure [T2l In figure [T3l 
the growth of the density perturbations for different values of the comoving wave 
numbers and k± is shown. In the direction of anisotropy the fc = mode initially 
decreases, but then turns and starts to grow unboundedly. The higher modes initially 
have decreasing amplitudes. Also these starts to increase, but more slowly, at about 
the same time as the fc = mode. In the perpendicular directions the modes do not 
show this initial decrease, but otherwise behave in a similar way. 

6.3.4- +K- to^dS, dust This dust background starts close to an expanding Kasner, 
and ends at an expanding de Sitter, -^dS. For the evolutions of the background 
quantities see figure 1121 for the corresponding radiation case, which shows a similar 
behaviour. In figure [H] the growth of density perturbations for fc = 1 is given. 

7. Summary 

A closed system for scalar perturbations on Kantowski-Sachs cosmologies has been 
found in terms of gauge invariant variables. For long wavelengths and low densities 
some analytical results are obtained. By redefining some apparently singular terms, 
the system could be rewritten in a form suitable for numerical analysis. Due to the 
complexity of the governing equations, the choice of background, initial conditions and 
wave numbers, many different behaviors of the growth of the density perturbations can 
be obtained. In general the growth of density gradients is different in the direction 
of the anisotropy and in the perpendicular directions, which should infiuence the 
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formation of structures. 

More importantly, growth rates are also affected by a bounce, see, e.g., figure 
m Typically the density gradients experience maxima (or minima) slightly after 
the bounce in the anisotropic Kantowski-Sachs models. A comparison with density 
growths in bouncing closed FLRW models, see [Appendix E shows that the amplitudes 
of the perturbations are roughly constant initially and that the local extrema around 
the time of the bounce are absent or very small in the isotropic case. 

It was shown in [68] that a FLRW model with dust, radiation and cosmological 
constant cannot undergo a bounce, unless it took place after the formation of the 
quasars. In section 13.51 we use similar arguments in the Kantowski-Sachs case to 
show that a bounce in the directions perpendicular to the preferred direction cannot 
take place at all, whereas a bounce in the preferred direction cannot be excluded by 
these arguments. We also derive relations between present values of observational 
quantities including expansion and deceleration rates in the different directions. Since 
anistropies in the Hubble and decelerations parameters cannot be excluded by present 
observations, see e.g. [70] and [71], it would be of interest to get better estimates of 
these parameters. 

The present analysis can be extended to study tensor perturbations, including 
gravitational waves, and the coupling between scalar and tensor perturbations. In this 
way possible observational consequences of anisotropic bounces can be formulated in 
a way which makes it easy to compare with what happens in the standard case. 
The techniques developed here can also be used to develop a general framework 
for describing the evolution of cosmological perturbations in other non-standard 
cosmologies for example the Lemaitre-Tolman-Bondi models [T^. Such studies are 
extremely important as they test the various rigidities which exist between different 
set of independent observables ( from the background and at the level of perturbation 
theory) that can be used to test the underlying hypothesis of the standard cosmological 
model. 
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Appendix A. Some approximate background-solutions with matter 



Here some approximate solutions with /i ^ p <C A to the system (|53p -(|55 p are given 
With the bounce solution ([49]) as background the solution is given by 

= Ccosh"''(\/At) 



Ml 
Si 



— cosh(vAt) X 



dt 



+ = cosh-2(VAt) X [E+ 

6 4 



cosh''+\^/At) 

cosh^-"'{VAt)dt 



(A.1) 
(A.2) 

A.3) 



where C, D and E are constants of integration. The integrals in (|A.2I) and (IA.3I) can 
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be performed for dust, 7 = 1, giving 

C C 

—= ta,-nh{VXt) and — -=smh{VAt) (A.4) 

VA 4vA 

respectively, and also for stiff matter, 7 = 2, giving 

(e^*)) and ^ (A.5) 



f sinh(VAi) /.yxA^ Ct 



VA y2cosli^(VAi) 
respectively. 

With the critical points ±A, (|45p . as background, the solutions are given by 

^^l = Ce^^^* (A.6) 

El - = Z?e±^* T — ^^e^^-^* (A.7) 
3 (7 + l)VA 

El + = i?eT2VAt ± (6 - 57)C ^^^VAt ^ (^.8) 
6 4(2-7)VA 

Note that since for all of these solutions some terms grow unboundedly, they are only 
valid for limited time intervals. 

Appendix B. Calculation of some first order quantities 

In the calculations several contractions of the type XabU""^ of PSTF tensors 

Xab = XiUaUb ~ ^Nab) + ^Xi^^nb) + Xab (B.l) 

will appear. To first order contractions of two, three or four PSTF tensor are given 

by 

XabV"'' ^\XY , XabV^'^Z^ ^-^XYZ , Xaby'^Z.dU'''' ^^XYZU (B.2) 

respectively, so that, e.g., = ^(Tabcr"^ = 

In calculating the expression for S = (o'^'^Sab)' , needed for the equation for Sa, 
the following term 

a"''(curliJ),b = (7"''(curl(curla)),, = a''''if''<a^ ,{cmla)b>d = 

a''''v"'a^c{cm\a)bd = I'^'v^'^a^crf^Kb^eTdyf = <J%V'daV'^<'^ e'T^'^f , (B.3) 

where it was used that aab already is PSTF and that the projected derivative of rjabc 
is zero, will appear. Using the definition of PSTF and that 





ud ue 
"a "-a 


"-a 




VabcV'''^ = 


hi hi 


h^ 

"■b 


(B.4) 




hi hi 


H 





it simplifies to 

= -V,V'^(a2) + la^b^bi^^aO) + aab^c^^a"' , (B.5) 

where was used in the last step. With the definitions ([7T|) and the 3-dimensional 
Ricci identity 

VcV^CTafa — ^d^cf^ab ebcdCl , (B.6) 
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that holds with Vq since ujab — 0, one gets 
a a 

The 3-Riemann tensor is given in terms of the projected 4-curvature, Rabcd± = 
^a^b^c^d^f'fghj s-iid the extrinsic curvature, Kab = '^{aUb) = 6*73 + aab, as 



^'^^Rabcd — Rabcd± 

and the 3-Ricci tensor by 



KacKbd + KbcKa 



^'^i?afc = Sab + ^hab ^'^i? = Sab + ^hab ( 2fl - -6 + 2a^ + 2A I .(B.9) 



Then, on using that the 4-curvature can be decomposed as 

R'^'cd = ^(m + 3p - 2A)uKi,h% + h%h% + 

4uKi,E% + 4h%E%] + 2i^'^'-u^,Ha]e + 2i^,deU^'^ H^^^^ (B.IO) 



(|B.7|) becomes 



1 



1^ 



a a 
where S is defined from equation (j74l) 



a''^{cm\H)ab = -o^^VaZb ~ -V aT" + 2 {(7"^) - -O^a^ + iSa^ + 2a'^{n + A) (B.f 1) 



Appendix C. Perturbative equations in orthogonal directions 

The perturbative equations, obtained by projecting ([76| -([79 |) with Nab^ are 



/i 



2 /i 



Z-a = (Za)- = -hz-a - ^flV-a + ]-^Z-a - 2% + ^ 

3 2 2 2 11 + p 



S+^T.^]Va 



Ta = {%)■ = -29% + -Sa~ ^Y.'^Z-a + ( 5 + 

fi + p 2 2 



^J. + P 



■Vr, 



(C.l) 



(C.2) 



(C.3) 



s-a-{s-ay^U' + 2l-]%+r-j:-le-2l;]s, 



E2 

5 X 



3 S 
p'fi 



M + P 



3 X 



f 



S - -OT. + 2S^ (JaVPb 



(C.4) 



The system obtained by taking the 2-divergence of (|C.ip - ()C.4p is 



/i 3 



2^ 2/i + pV 2 



PP. ^2 



+p 



P + ^ 



(C.5) 



(C.6) 
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S M + P 



-I- 



fip' 1 



Ai + P 



25^ 7" 



£3 
S 
lip 



M + P 

Appendix D. Perturbations around analytical solutions 

Appendix D.l. Perturbations around -^X 

With the +X-solution as background the system (P51) - ([M)) for the hat- variables reduces 
to 



(C.7) 



(C.8) 



I? = - -\/A2? - Z , Z = --\/liZ - 2f 
3 3 

3 3 ' 27 



22 



-AT 



(D.l) 
(D.2) 



if the harmonic numbers are put to zero (corresponding to large wave-lenghts of the 
perturbations). From these equations one finds the fourth order equation (|109p for the 
density perturbations T) with solution (jllOp . The other quantities are then obtained 
as 

35 



givmg 



i = - 2? - , f =-'b + 2VA + —KV 

3 ' 2 18 

5 = _ ip(3) _ ^VAf) - — - ^A3/2p 
2 6 18 54 



'b = Ai+A2 + Ai, i = VA(2Ai - A2) - A3 



T 



A 



(4Ai 



VA 



A, 



5 



A3/2 



(4Ai 



11 



A2) 



as intial conditions at t = 0. 



(D.3) 
(D.4) 

(D.5) 
(D.6) 



Appendix D.2. Perturbations around vacuum bounce solution 



With the vacuum bounce solution pQ]) as background the system (|93l) -(|96 |) for the 
hat- variables reduces to 

5, 



V 

z 



-ev- z 
3 

I.i-Vt + 2^5 



(D.7) 
(D.8) 
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i ^ (5A - 'f') i + ^ (a= + le-) f - 2« (l + i (D.IO) 

(once again only considering the long wave-length limit). It is suitable to change the 
independent variable to 9 through 

The equation (|112p for V with solution (|113l) is then obtained. The other quantities 
are given by 

-lev-iA-e')^ (D.ii) 

'^^Wa{1 (9A2 - _ 8^4) p_ 5^ ^ 4^2) - 0^) ^ + 

3 (3A - 56^) (A - ^ - 96 (A - 0^)^ (D.12) 
S =-Ap0(57A + 802)_ JL^ (33 A + 20 (A _ _ 

Taylor expanding 2? to second order around 9 — gives 

t)^Ai + {A2+A3)9-^^9^ + ... (D.14) 
6 A 

i = - (^2 + i3)A + (i^i2 + ^is)^?' + ... (D.15) 
6 3 

t^-fi.-?i3. + ^^^^ + ... (D.16) 
5 = (^^2 - Ia3)A' - 3A^A9 - (^i2 + ^i3)A02 + (D.17) 
where Ai = AiA^/e, A2 = v42A^/6, ig AaA^s and i4 = A4A5/6. 

Appendix E. Perturbations on bouncing closed Priedmann 

Closed Friedmann models with a positive cosmological constant may undergo a bounce 
The evolution equations are given by 

fi=-9{fi + p) (E.l) 

^ = - J^'-^(A* + 3p-2A). (E.2) 

With an equation of state p = (7 — the following integral 

6|2 = S^i + 3A - 3Cn^/^'' (E.3) 



is found. Here closed models correspond to C > 0. In figure [Ell a bouncing radiation 
solution is shown. It starts from a contracting state, goes through a bounce and then 
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Figure El. A closed initially contracting Friedmann model with radiation, tliat 
goes through a bounce and then starts expanding. Initial values at to = 1 are 
given by fio = 0.001, Bo = -1.65. A = 1 and C = 2.96. 



asymptotically approaches de Sitter. 

To study the growth of density perturbations it is sufficient to use the density 
gradient T^a = "'^"^ and expansion gradient Za — aS/ a(^, in which the system now 
closes [ST]: 

2?a = -I?a-2,(l + ^) (E.4) 

t^ 

Za = -IdZa - i/iPa + f-^ - A + - a^'V, . (E.5) 

3 2 + 6 J M + P 

Here p' = dp/dfi. In accordance with equations (pi)) and we define the scalar 
density gradient as 

V = aV^Pa • (E.6) 
A harmonic decomposition 

T^ = J2'^kQk wliere V^g^ = -Qj^ and (9^ = (E.7) 

k 

is then done. Here fc^ = n{n + 2) for n = 1,2, [7^175] . 

The growth of the density perturbations for the radiation background in picture 
lEll is shown for different wave-numbers in picture IE2I Initially the amplitudes of 
the perturbations are roughly constant with a small dip in the amplitude of the 
fc = mode, that serves as an approximation for the long wavelength limit, before the 
bounce. 
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Figure E2. The growth of the density perturbations D in the background given 
by figure [ET] for the wave numbers k = k/ag = 0, 1, 5 and 20. Initially, at io = 1, 
V = 0.01. The time of the bounce is indicated with a dotted vertical line. 




Figure E3. The growth of the density perturbations V for a dust background. 
Initially, at to = I, "D = 0.01. Initial conditions for background fio = 0.001, 
0Q = —1.65. A = 1 and C = 9.35. The time of the bounce is indicated with a 
dotted vertical line. 
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In the corresponding dust case the behaviour of the background quantities are 
similar to the radiation case, whereas the evolution of the perturbations is independent 
of the wave- number, as can be seen from (jE.5|) . See IE3I for the growth of V . 

Bouncing closed Friedmann universes containing dust and radiation with 
cosmological constant are excluded by observations, due to an argument by Borner and 
Ehlers '68] . We here repeat their argument for the dust case in terms of our variables. 
Using the definition of the scale factor a in terms of the expansion, 9 = Sa/a, and Eq. 
(|E.1|1 one obtains /i = /io(ao/a)'^, where the subscript refers to present values. The 
constant in (jE.3[) is given by 3C/iQ^'^ — 3/io + 3A — in terms of present values. At 
the bounce we have d* = and d* > 0. Using this into Eqs. (|E.2p and (jE.3[) one then 
obtains the following inequality 

f7„ [(z + l)3-3(z + l) + 2] < 2 (E.8) 

where = 3/io/^o ^^'^ present fraction of dust relative to the critical density and 
z = ao/ai, — 1 is the redshift corresponding to the time of the bounce. If the bounce 
took place before the formation of quasars, that are observed at redshifts z > 4, flm 
cannot overcome the value 0.02, in contradiction with observations. 
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